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■ We study the nuclear liquid-gas phase transition on the basis of a two-component 

lattice gas model. A Metropolis type of sampling method is used to generate micro- 



scopic states in the canonical ensemble. The effective equation of state and fragment 
mass distributions are evaluated in a wide range of temperatures and densities. A 
definition of the phase coexistence region appropriate for mesoscopic systems is pro- 
^ ' posed. The caloric curve resulting from different types of freeze-out conditions are 



presented. 



It is commonly believed that nuclear matter undergoes a liquid-gas phase transition 
at lower densities and temperatures. In the early 1980s it was suggested |l|] that this 
phase transition might be probed in heavy-ion collisions at intermediate energies by 
observing the disintegration of colliding nuclei into many fragments of different sizes, a 
phenomenon commonly known as multifragmentation. 

A wide variety of models has been proposed for nuclear multifragmentation, ranging 
from statistical to dynamical ones (see e.g. Q for a review). In principle, most of these 
models yield an equation of state of nuclear matter. The evaluation of the equation 



of state is, however, in general rather difficult and has only been done in some simple 
cases [[|. In this paper we focus on the relationship between the fragmentation of 
a nuclear system and the corresponding equation of state, calculated on the basis of 
a two-component lattice gas model. In particular, we will discuss the nuclear caloric 
curve (temperature versus excitation energy) resulting from different types of freeze- 
out conditions. The details of the model and the numerical simulations will be further 
discussed in a forthcoming article 0]. Here, we only present some main features and 
results of the model. 

Common for all types of lattice gas models is the use of a spatial lattice, with a 
certain dimension, number of sites, Ny, and a certain Bravais-structure, such that each 
lattice point has a fixed number of neighbours 7. For simplicity and for comparison 
with earlier works we consider here a simple cubic (s.c.) lattice with 7 = 6. The lattice 
spacing a$ is chosen such that the density of a fully occupied lattice equals the nuclear 
ground state density po = 0.16 fm~ 3 , i.e. a[j = 1/po- The effective nucleon-nucleon 
(iV./V)-interaction is represented by a nearest-neighbour square well potential V with a 
repulsive hard core. In terms of the discritized lattice distance r between the nucleons 
(r = 0, do, \/2ao, • • •) the potential is given as, 
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00 for r = 
e t3 . t3j for r = a 
for r > ao 



Here, £3. denotes isospin of nucleon i and et 3i t 3j are the interaction strength parameters. 
In fact we only allow for two such parameters, e s for pp and nn interacting pairs and 
€d for pn pairs. According to the dynamical potentials and previous works on lattice 
models || is chosen slightly repulsive, = 1 MeV. The other parameter, e s , is chosen 
such that the nuclear binding energies are reproduced, e s = —5 MeV. 

If p~i, fl denote respectively the momentum and position of nucleon i, the Hamiltonian 
of the system of A nucleons can be written as, 

w = r + v = Ef£ + £W). (i) 

i i<j 
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where m is the nucleon mass. It can be shown that the interaction Hamiltonian V 
is homomorphic to the spin-1 Ising model, which is well-known in solid state physics. 
If e s is put equal to the model reduces to the standard lattice gas model which is 
isomorphic to the spin-1/2 Ising model. This model was first considered in context of 
nuclear fragmentation in || and later used by S. Das Gupta and J. Pan []|, ||, ^] and X. 
Campi and H. Krivine J|]. The two-component lattice gas model has also been used in 

SI- 

In the following we will for simplicity treat the model in a canonical ensemble. This 
implies, that the temperature T, the volume V, the number of particles A and the 
asymmetry of the system, y = Z/A, are the control variables. In this paper we only 
consider symmetric matter, y = 0.5. 

The states / of the nuclear system can be classified by the positions, momenta and 
isospins of all the nucleons, / = {(r, p, £3)4 |i = Due to the classical nature of 

the model the canonical partition function, Z, can be factorized into two parts coming 
respectively from the momentum space and the configuration space. The first part is 
proportional to the partition function of an ideal gas. The contribution of the momenta 
of the nucleons to the various thermodynamical quantities can therefore be calculated 
in a straightforward way. 

The partition function and thus the Helmholtz free energy, F, of the configuration 
space is far less trivial to obtain. In this work it is done numerically by performing a 
Monte Carlo type of sampling. The realization of this sampling method consists of three 
parts: the initialization, the equilibration of the system and the calculation of ensemble 
averages. 

The initialization of the system is done by putting the particles one by one in the 
lattice. At each step the spatial probability distribution is changed with a Boltzmann 
factor, e~ ei il T , in the neighbourhood of the first particle. This method of initializing 
the system is identical to the one used in ||. 

The equilibration of the system is performed by constructing a random walk of 
microscopic states {f u } through the configuration space via a Markov proces^, starting 

A sequence of states each of which depends only on the preceding one 
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from the initial state, f u =o- The Markov process is defined by specifying a general 
transition probability matrix, p(f — * /') from the state / to the state /'. The choice of 
the transition matrix is dictated by the detailed balance condition. In the present work 
it is defined as, 

-/9AV ifA v>0 

(2) 

1 otherwise 

Here, the test function, /', is obtained from / either by changing the position of a ran- 
domly chosen nucleon or by interchanging a randomly chosen proton-neutron pair. In 
case of acceptance, the trial function /' will be the new state of the system f u +i = /'• 
Otherwise, the system remains in the same state, f u +i = fv It can be proven for- 
mally |l(J that the above choice of p leads to the desired convergence property, i.e. the 
equilibrated probability distribution of microscopic states will correspond to a macro- 
scopic state of minimum free energy. This algorithm is similar to the one introduced by 
Metropolis et. al £□}. 

When equilibrium has been reached, the ensemble average < Q > of any quantity 
Q can be calculated by repeating the same procedure a certain number of times and 
recording the value of the quantity, Q u , at each time step v. 

By a systematic variation of the density p = A/V and the temperature T the algo- 
rithm one can eventually probe all the different points of the phase diagram including 
the coexistence region. The density is varied by varying the number of particles in 
the lattice for a fixed volume. The volume is put to Ny = 6 3 which approximately 
corresponds to the maximum fragment mass, A max = 190, observed in heavy ion ex- 
periments. The boundaries of this lattice are chosen to be open to mimic the effect 
of the vacuum surrounding the real nuclear system. The results for periodic boundary 
conditions will be presented elsewhere 

The quantities, Q, of interest are the ensemble average E(T,p) =< V > of the 

potential energy V, and the distributions of cluster sizes, i.e. the mean number of 

fragments, or the yield Y(A) as a function of the mass number A. Contrary to most 
2 This is an important difference to the grand canonical ensemble which only allows the system to be 
in a single phase. 



P(f - /') = 
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statistical models of nuclear fragmentation (see i.e. 0] and references herein), clusters 
are not explicitly denned in the lattice gas model. In our simulations we implement a 
definition based on pairwise binding. Two neighbouring nucleons are thus prescribed to 
the same cluster if their relative kinetic energy is not sufficient to overcome the attractive 
bond: <j 2 /4m + et 3i t 3 . < 0, where q = P\~P2 is the relative momentum between the two 



particles. This definition was originally proposed by Hill \12\ and has also been used by 
S. Das Gupta et. al ]|]. 

By applying the standard laws of thermodynamics all other quantities can be de- 
rived from the internal energy E of the system in the different points, p, T, of the phase 
diagram. One of the most interesting applications is the study of the liquid-gas phase 
transition in a finite system. The first order phase transition is signalled by the presence 
of the phase co-existence region. In a finite system the definition of this region is some- 
what ambiguous. In the calculations below such regions are identified by the existence 
of two distinct densities, p g < pi, for which the following conditions are fulfilled, 

pi(pi,T) =p g ( Pgi T) = p, F{p h T)/V = »{p l -p g ) + F{p g ,T)/V (3) 

The points (p g ,T) and (pi,T) belong to respectively the gas branch and the liquid 
branch of the co-existence curve. Geometrically, the above criteria correspond to the 
requirement that the points (p g , f g ) and (pi, fi) belong to a common tangent of the free 
energy density, / = F/V . For a macroscopic system the free energy density would 
follow a straight line in the region p g < p < p\. For a mesoscopic system the free energy 
density is slightly larger in this region due to the presence of an interface between the 
two phases. 

The equation of state for the system with volume Ny = 6 3 is depicted in Fig. 1 with 
the pressure P as a function of density p/po and temperature T. The thick solid line is 
the co-existence curve calculated according to the definition (|3|). The two branches of 
the curve terminate at the critical point (p/ po,T,p) cr ~ (0.39, 4.8 MeV, 0.15 MeV/fm 3 ). 
For comparison, the critical point in the macroscopic limit is found to be (p/po, T, P) cr ~ 
(0.58,6.7 MeV, 0.36 MeV/fm 3 ) ]J. In the co-existence region for a fixed temperature 
the pressure as a function of the density is almost a constant. 
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One should bear in mind that in experiment one observes only a final state of the 
system where the fragments cease to interact (freeze-out configuration). It is not clear 
a priori which thermodynamical conditions correspond to these configurations. Below 
we consider two possibilities corresponding to constant pressure and constant density at 
freeze-out. The two lines across the co-existence region are defined by p = 0.10 MeV/fm 3 
and p/po = 0.3, respectively. The letters A — E along the isobar are chosen with a fixed 
step in densities Ap/po = —0.15 starting from p/po = 0.75. the letters along the isochor 
are chosen at the temperatures T = (2.4, 3.5, 4.3, 4.9, 6.5) MeV (points c and D coincide). 

The cluster size distribution corresponding to these points of the phase diagram are 
shown in fig. 2a and 2b for the isobar p = 0.10 MeV/fm 3 and the isochor p/po = 
0.3 respectively. At point A the system essentially consists of a compound nucleus 
surrounded by nucleons and small clusters. Point A belongs to the borderline of the 
co-existence region. At point B the compound nucleus starts to break-up producing 
some intermediate mass fragments (IMFs: 4 < A < 20). At point C and D fragments of 
all different sizes are present, whereas at point E the size distribution is exponentially 
decreasing indicating that only the gas phase is present. Thus the appearance of IMFs 
is a clear manifestation of the first-order liquid-gas phase transition. 

The evolution of the size distribution is observed to be smoother along the isochor. 
Point a belongs to a transition region for which the pressure increases rather abruptly. 
At even lower temperatures the system essentially consists of one compound nucleus 
surrounded by vacuum (not shown). The sudden increase of the pressure results from 
the break-up of this nucleus. As can be seen in the plot this produces some IMFs. The 
distribution at point d still contains a significant contribution of IMFs although it does 
not belong to the co-existence region, according to definition of eq. (jH). Thus, the cluster 
size distribution alone can not give a unique relation between multifragmentation and 
liquid-gas phase transition. At point e the distribution is closer to the exponential one. 

Figures 2a and 2b indicate that cluster size distributions are widest at the vicinity 
of the critical point. If the cluster size distribution is approximated by a power law 
Y(A) ~ A~ T , one finds that a good fit can be achieved for a wide range of temperatures 
and densities around the critical point. The range of fragment masses used for the power 
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law fit was chosen to 1 < A < 20, but the conclusions are not very sensitive to this choice. 
Typical values of x 2 P er degree of freedom are in the range of 0.03 — 0.13. Our analysis 
shows that for each isochor, r has a minimum at the binodal line. Furthermore, r has 
a global minimum point in the vicinity of the critical point, where r ~ 2.2. This is in 
accordance with the value of the Fisher exponent for the Ising model, r = 2.21. 

From the yield of clusters Y(A) it is not easy to distinguish between the two types of 
freeze-out conditions. The difference reveals itself more clearly in the associated caloric 
curves, showing the variation of the temperature with respect to the excitation energy 
E* per nucleon. Here, the excitation energy is defined as 

E*(A,T) = 3/2 -T ■ (<M> -1) + E(A,T) - E(A,T = 0), 

where <M>= J2aY(A) is the total multiplicity of the system. Note, that we identify 
the kinetic part of the excitation energy with the total kinetic energies of the clusters 
rather than the total kinetic energy T of the individual constituents. The caloric curves 
for p = 0.1 MeV/fm 3 and for p/po = 0.3 are shown in Fig. 3. It is seen that the 
pronounced plateau in the caloric curve appears only in case of the fixed pressure con- 
dition. In both cases the caloric curve approaches a straight line for high temperatures. 
At low temperatures the model predictions are not accurate, since quantal effects cannot 
be neglected in this region. Qualitatively, the caloric curve at constant pressure looks 
very similar to the one predicted by the Statistical Multifragmentation Model 13] and 



recently observed in the ALADIN experiment [14]. 

To each of curves on the phase diagram, which correspond to a specific change of 
the system state in a reversible process, a particular heat capacity C can be attributed. 
Our analysis shows, that the specific heat capacity for fixed pressure, C p , in fact has 
a singularity along subcritical isobars. This infinity is related to the latent heat of the 
transition. On the other hand, due to finite size effects the specific heat capacity for 
fixed number of particles and volume, Cy, is finite in all points of the phase diagram. 
We will elaborate further on these points in a forthcoming publication 

In conclusion, our results demonstrate the close connection between the nuclear 
liquid-gas phase transition and the nuclear multifragmentation. The calculations have 
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been made on the basis of a two-component lattice gas model. A sampling method has 
been presented by which both the effective equation of state and the fragment mass 
distributions have been evaluated in a consistent way. The method has an advantage 
compared to other statistical approaches to nuclear multifragmentation, where the link 
between the cluster observables and the corresponding thermodynamical properties of 
the nuclear system is less transparent. 

We have proposed a definition of the phase transition, more relevant for finite sys- 
tems, which is based on the identification of an interphase between the two phases. The 
relevance of this definition is supported by the critical exponent analysis and the singular 
behaviour of the specific heat capacity, C p , along subcritical isobars. We have examined 
the caloric curve resulting from two different freeze-out conditions. It was demonstrated 
that the occurrence of a plateau in the caloric curve seems to be consistent with the 
condition of constant and subcritical pressure in the freeze-out configurations. 

There are of course some obvious deficiencies of the model. First, the true Fermi- 
statistical nature of the nucleons can not be neglected for lower excitation energies. 
Secondly, the idea of a thermal bath is certainly an idealization and it is necessary to 
perform a real microcanonical sampling in the Monte Carlo algorithm. Third, for larger 
systems it is necessary to include the Coulomb repulsion and determine the correspond- 
ing change of the phase diagram. Finally, the prescription of clusters in the model is still 
not unambiguous. The clusters defined by the pairwise binding are generally excited 
and should undergo de-excitation at later stages (see (ill). In the future, we plan to 
improve the model in various respects to take these deficiencies into account. 
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Figure captions 

Fig. 1. Pressure P as a function of density p/po and temperature T for a lattice with 
Ny = 6 3 sites. The thick solid line is the co-existence curve. The line with capital letters 
is the isobar p = 0.1 MeV/fm 3 . The line with small letters is the isochor p/po = 0.3. 

Fig. 2. The yields Y(A) as functions of the fragment mass A corresponding to 
different points of the phase diagram: a) The points (A — E) belong to the isobar 
p = 0.10 MeV/fm 3 , b) the points (a — e) belong to the isochor p/po = 0.3. Note the 
different scales on the x-axis. 

Fig. 3. The temperature T as a function of the specific excitation energy E*/A for 
two different freeze-out conditions. The solid line corresponds to the fixed pressure 
condition, p = 0.1 MeV/fm 3 , and the dotted line corresponds to the fixed density 
condition, p/po = 0.3. The letters along the lines refer to fig. 1. 
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Figure 1 
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Figure 2a 




Figure 2b 
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Figure 3 
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